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Abstract We analyzed the space velocities of Gould Belt stars younger than 125 Myr 
located at heliocentric distances < 650 pc. We determined the rotation and expansion pa¬ 
rameters of the Gould Belt by assuming the existence of a single kinematic center whose 
direction was found to be the following: lo = 128° and Ro = 150 pc. The linear velocities 
reach their maximum at a distance of ~ 300 pc from the center and are —6 km s“^ for 
the rotation (whose direction coincides with the Galactic rotation) and -|-4 km s“^ for the 
expansion. The stellar rotation model used here is shown to give a more faithfid description 
of the observed velocity field than the linear model based on the Oort constants Aq and 
Be- We present evidence that the young clusters [3 Pic, Tuc/HorA, and TWA belong to the 
Gould Belt structure. 

Key words: Gould Belt, Galaxy (Milky Way), OB associations, TW Hydrae, [3 Pictoris, 
Tucana. 


INTRODUCTION 

The star-and-gas complexes associated with star formation processes in our Galaxy and 
in other galaxies have a hierarchical structure (Efremov 1998; Efremov and Elmegreen 1998). 
Giant star-and-gas (GSG) complexes trace the spiral pattern of the Galaxy (Efremov 1998). 
The GSG complexes have masses of ~ 1 x lO^Mo, sizes as large as 1000 pc, and lifetimes 
r < 10® yr. Less massive structures, such as OB associations, open star clusters, and giant 
molecular clouds, are part of a GSG complex. The Sun lies within a GSG complex that is 
known as the Gould Belt (with a radius of ~ 500 pc and a lifetime of r 60 x 10® yr). In 
turn, the Gould Belt is part of an older (r ~ 5 x 10® yr) and more massive 2 x lO^M©) 
structure about 1000 pc in size that is known as the Local (Orion) Arm or the Local system 
of stars perceived as a gravitationally bound, long-lived system. Its kinematics suggests that 
it is associated (Olano 2001) with the Sirius star supercluster (Eggen 1984, 1992). Since the 
mean residual velocity of the Sirius supercluster stars relative to the local standard of rest 
(LSR) is low, the center of mass of the Local system of stars has made several turns around 
the Galactic center in a nearly circular orbit. Over its lifetime, the Gould Belt has probably 
experienced a single impact from a spiral density wave. 

An important indicator of the stability of such a system as the Gould Belt is the existence 
of proper rotation. Such rotation is reliably found from observations of Gould Belt stars 
(Lindblad 2000; Bobylev 2004). Over its lifetime, the Gould Belt has been affected by the 
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shock wave generated by a spiral density wave at least once. An analysis of the motions 
of the OB associations that fill the interarm space (r < 3 kpc) reveals a complex periodic 
strnctnre of their residnal velocities attribntable to the inflnence of spiral density waves 
(Mel’nik et al. 2001; Zabolotskikh et al. 2002; Mel’nik 2003). One of these featnres is the 
prominent structure in the radial residual velocity distribution in Galactocentric distance 
in the Gould Belt region (see Fig. 3a in Mel’nik et al. 2001) that we associate with the 
well-known positive K-effect or with the expansion of young stars. In this paper, based on 
a sample of individual stars that are members of OB associations, we study this effect in 
detail. 

The goal of this work is to study the motion of the Gould Belt using data on stars with 
reliable age estimates. For this purpose, we use young and nearby OB associations with 
age estimates obtained by different authors on the basis of currently available observations 
and methods. The rotation and expansion parameters of the Gould Belt are determined 
by assuming the existence of a single kinematic center. To this end, we use Bottlinger’s 
formulas in a form that allows the direction of the single kinematic center to be analytically 
determined simultaneously with the rotation and expansion-contraction parameters. 

METHODS OF ANALYSIS 

Bottlinger’s Formulas 

In this paper, we use a rectangular Galactic coordinate system with the axes directed 
away from the observer toward the Galactic center (/ = 0°,6 = 0°, the X axis), along 
the Galactic rotation (/ = 90°, 6 = 0°, the Y axis), and toward the North Galactic Pole 
(6 = 90°, the Z axis). We derived the basic equations from Bottlinger’s standard formulas 
(Ogorodnikov 1965). By assuming the existence of a single kinematic center of rotation and 
expansion- contraction (G, i?o), we transformed the formulas to 

Vr = ug cos b cos I + vg cos 6 sin / -|- wg sin 6 — (1) 

—Di{R — Ro) sin(/ — G) cosb— 

—D 2 {R — Ro) cos(/ — lo) cosb— 

—Fi {R — Ro)^ sin(/ — /„) cos b— 

—F 2 {R — Ro)^ cos(/ — lo) cos6-1- 
YkoT cos^ b + k'^r{R — Ro) cos^ b+ 

+0.5k'^r{R — Ro)^ cos^ b, 

4:.74:rfj,i cos b = —ug sin I + vg cos 1+ 

+uJor cos b + uj'j-{R — Ro) cos b+ 

+t)Roj'lr cos b{R — Ro)^— 

—Di{R — Ro) cos(/ — lo) + D 2 {R — Ro) sin(/ — lo) — 

-Fi{R - Ro)'^ cos(/ - G) + F 2 {R - i?o)^sin(/ - G), 

4.74r/rb = —ug cos / sin 6 — vg sin I sin b + wg cos b+ 

+Di{R — Ro) sin(/ — lo) sin6-|- 
+D 2 {R — Ro) cos(/ — lo) sin6-|- 
+Fi{R — /?o)^sin(/ — lo) sin6-1- 
+F 2 {R — RoY cos(/ — lo) sin6— 
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—koT cos 6 sin 6 — k'j-{R — Ro) cos b sin b— 

—0.5k”r{R — RoY cos & sin b. 

Here, the coefficient 4.74 is the quotient of the number of kilometers in an astronomical unit 
by the number of seconds in a tropical year, r = I/tt is the heliocentric distance of the star, 
Ro is the distance from the Sun to the kinematic center, and R is the distance from the 
star to the kinematic center. Since we corrected the observed motions for the peculiar solar 
motion with respect to the LSR (Dehnen and Binney 1998), the quantities uq, vq, and wg 
are (in contrast to the standard approach, we reversed the signs) the velocity components of 
the centroid of the stars under consideration with respect to the LSR; the components of the 
stellar proper motion nicosb and fib are in mas yr“^ (milliarcseconds per year), the radial 
velocity R- is in km s“\ the parallax tt is in mas (milliarcseconds), and the distances R, Ro, 
and r are in kpc. The quantity Uo is the angular velocity of the stellar system at distance Ro; 
ko is the radial expansion (or contraction) rate of the stellar system at distance Ro; ui”, 
and kY /c" are the corresponding derivatives with respect to the heliocentric distance taken 
at distance Ro; and lo is the direction of the kinematic center. R can be calculated using the 
expression 

R^ = (r cos bY — 2Ror cos b cos(/ — k) + Rl- (2) 

In Eqs. (1), we use the unknowns 

Di = Dcoslo, D2 = —Dsml0, (3) 

Fi = F cos Ig, F2 = —Fsinle, 

where D = Df + D2 = ui'Ro m F = Ff + F2 = 0.5cc;"i?o- Based on relations (3), we 
determine the correction Ig to the assumed value of lo twice: 

tail{lg)D = —D 2 ID 1 , (4) 

tcm{lg)F = -F2/F1, 

in this case, the new direction is lo + le- On the other hand, we may introduce similar 
quantities, G = k'Ro and F[ = Q.bF'Ro, with the unknowns Gi, G2, Hi, and H2 for the 
expansion to obtain 

Gi = GcosIr = —G sin Ig, 

G2 = —G sin Ir = G cos Ig, (5) 

Hi = H cosIr = —H sin Ig, 

H2 = —HsinlR = H cos Ig. 

The orthogonality of the expansion and rotation effects {Ig = Ir + 90°) is reflected in relations 
(5); as a result, the unknowns in Eqs. (1) cannot be separated. Therefore, to determine 
the direction of the kinematic center, it will suffice to use one of the sets of unknowns, 

for example, Di, D2, Fi, and F2, in Eqs. (1). Based on the described approach, we can 

independently estimate the distance from the Sun to the kinematic center. This estimate 
can be obtained from the formulas that follow from relations (3)-(5): 

D = Ro\u'^ - k'J, ( 6 ) 

F = 0 . 5 Ro\u;”-k'Y. 
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Equations (1) contain thirteen sought unknowns: u©, Vq, Wq, Uo, K, k'^, k”, Di, 

D 2 , Fi, F 2 , which can be determined by the least-squares method. If the direction of the 
kinematic center is known, then the eqnations have the original form 

Vr = Ug cos b cos I + (7) 

+vg cos 6sin / -|- vjg sin b — Ro{R — Ro) x 
X sin(/ — lo) cos 6a;' — 0.5i?o(.R — Ro)^x 
X sin(/ — lo) cos 6a;'' -|- cos^ bkor+ 

+ {R — Ro){rcosb — Ro cos(/ — /o))x 
X cos bk'o + 0.5(i? — Ro)‘^{r cos 6— 

—Ro cos(/ — lo)) cosbk'o, 

4.74r/i/cos6 = —ucsin/-!- 
+vg cos I — {R — Ro){Ro cos(/ — lo) — 

—r cosb)uJo — 0.5{R — Ro)^{Ro cos(/ — lo) — 

—r cos b)Uo + r cos buo + Ro{R — Ro)x 
X sin(/ — lo))ko + 0.5Ro{R — Ro)^x 
X sin(/ - lo))ko, 


4.74r/ife = —Ug cos I sin 6— 

—vg sin/ sin 6 -|- wg cos 6 + Ro{R — Ro)x 
X sin(/ — lo) sin 6 a;' -|- 0.5Ro{R — i?o)^x 
X sin(/ — lo) sin 6 a;'' — cos 6 sin bkoV- 
— {R — Ro){r cos b — Ro cos(/ — /o)) x 
X suibk'o — 0.5(7? — i?o)^(r cos 6 — 

— Ro cos(/ — lo)) sm.bko. 

These equations contain nine sought unknowns: u©, u©, m©, Uo, a;', ui'l, ko, k'o, k”. A pecu¬ 
liarity of this method is that it requires the existence of derivatives only with respect to the 
distance. A sample of stars nniformly distribnted even in one Galactic qnadrant can satisfy 
this reqnirement. 

The Statistical Method 

We use the standard statistical method (Triimpler and Weaver 1953; Parenago 1954; 
Ogorodnikov 1965) that consists in determining and analyzing the symmetric moment ten¬ 
sor or the stellar residnal velocity dispersion tensor. When both the radial velocities and 
proper motions of stars are nsed to determine the six nnknown components of the dispersion 
tensor, we have six equations that can be written for each star. The semiaxes of the residual 
velocity (Schwarzschild) ellipsoid that we denote by <Ji^ 2,3 can be determined by analyzing 
the eigenvalnes of the dispersion tensor. We denote the directions of the principal axes of 
this ellipsoid by /i 2,3 and 61 2 , 3 . A pecnliarity of the approach nsed here is that the stellar 
velocities corrected for the pecnliar solar motion with respect to the LSR and for the general 
Galactic rotation are used as the residual velocities. 
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The Ogorodnikov-Milne Model 


In the linear Ogorodnikov-Milne model, we nse the same notation that was introdnced by 
Clube (1972, 1973) and used by du Mont (1977, 1978). We modify the linear model developed 
to describe the general Galactic rotation to describe the peculiarities of the velocity field of 
nearby stars. To within terms of the first order of smallness r/Ro 1, the observed velocity 
V(r) of a star with a heliocentric radius vector r is described by the vector equation 

V(r) - Ygr - yQLSR = Vg + Mr + V. ( 8 ) 


Here, Vg/? is the systematic velocity of the star attributable to the general Galactic rotation, 
V lsr is the peculiar solar motion with respect to the LSR, Vg(mg, vg, wg) is the velocity of 
the centroid of the stars under consideration relative to the LSR, V' is the residual stellar 
velocity (the residual stellar velocities are assumed to be distributed randomly), and M is 
the displacement matrix that describes systematic motions similar to the proper rotation and 
expansion-contraction effects. The components of the matrix M are the partial derivatives 
of u(mi,M 2 ,M 3 ) with respect to r(ri,r 2 ,r 3 ): 

(P,g = 1,2,3). (9) 


The matrix M can be separated into symmetric, M+, and antisymmetric, M~, parts. Fol¬ 
lowing Ogorodnikov (1965), we call them the local deformation and local rotation tensors, 
respectively: 


M+ = 

pq 


M 


pq 


1 / dUg \ 

2[drg drjj 

1 / dUp dUq\ 
2 [drg 

{P,(l = 1,2,3). 


( 10 ) 


The basic equations are 


Vr = ug cos b cos I + vg cos 6 sin / -|- wg sin b + ( 11 ) 

-|-r(cos^ 6 cos^/Mil + cos^ 6 cos/sin/Mi 2 + 

+ cos b sin b cos /M 13 -|- cos^ b sin I cos /M 21 -I- 
-|- cos^ b sin^ IM 22 + cos b sin b sin IM 23 + 

+ sin b cos b cos /M 31 + cos b sin b sin /M 32 + 

-hsin^ &M33), 


4.74r/i; cos b = —ug sin I + vg cos /+ 
-|-r(— cos 6 cos / sin /Mu — cos/>sin^ /M 12 — 
— sin/^sin/Mi3 -f cosbcos^ IM21 + 

+ cos b sin / cos IM 22 + 

+ sin/>cos/M 23 ), 
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4.74r/ife = —ug cos / sin & — vq sin / sin 6 + wg cos 6 + 

+r(— sin 6 cos 6 cos^ IMu — 

— sin b cos b sin I cos IM 12 — 

— sin^ b cos /M13 — sin b cos b sin I cos IM 21 — 

— sin b cos b sin^ IM 22 — sin^ b sin IM 23 + 

+ cos^ b cos /M 31 + cos^ b sin /M 32 + 

+ sin 6 cos bM^s), 

where the stellar velocity components corrected for the general Galactic rotation and for the 
peculiar solar motion with respect to the LSR appear on the left-hand sides. Equations (11) 
contain twelve sought unknowns: the three velocity components Vg{ug, Vg, Wg) and the nine 
components Mpg that can be determined by the least-squares method. The deformation and 
rotation tensor components can be calculated using the values of Mpq derived from relations 
(10). This method requires the existence of all derivatives with respect to the coordinates. 
A sample of stars that densely fill each coordinate axis is needed to meet this requirement. 

OBSERVATIONAL MATERIAL 

Data on Stars 

We compiled our working list of stars located no farther than ~ 650 pc from the Sun 
using the following sources: we took the Hipparcos (ESA 1997) star numbers for the Col 121, 
Per OB2, Vel OB2, Tr 10, LOG, UCL, US, Cep OB2, Lac OBI, Cep OB 6 , f3 Per (Per OB3), 
and Cas Tau associations from de Zeeuw et al. (1999); for the Ori OBI a-d associations 
from Brown et al. (1994); for the IC 2391, IC 2602, NGC 2232, NGC 2451, NGC 2516, NGC 
2547, and Pleiades clusters from the list by Robichon et al. (1999); the list of members of 
the TW Hydrae (TWA), f3 Pic, and Tucana/HorA associations from Song et al. (2003); and 
for the a Car cluster from Platais et al. (1998). The data on 24 stars having X-ray emissions 
and located no farther than 50 pc from the Sun were taken from Makarov (2003), where 
they are denoted as an XY sample. The star HIP 30030 in Makarov’s list is a member of 
the TWA cluster, belongs to the Tucana/HorA group, and is designated TWA42 (Song et 
al. 2003). The data on 34 stars located no farther than 100 pc from the Sun were taken from 
Wichmann et al. (2000, 2003). Pour stars were not included in other lists; we placed them in 
group 2. Table 1 gives the ages of the selected clusters. We took the equatorial coordinates, 
parallaxes, and proper motions from the Hipparcos catalog and the radial velocities from 
the catalog by Barbier-Brossat and Pigon (2000). We use only single stars (the astrometric 
orbital binaries marked by the symbol O were rejected) for which the parallaxes, radial 
velocities, and proper motions are available. Based on the stellar ages (Table 1), we divided 
the stars into three groups: (1) the youngest stars with ages < 10 Myr (group 1); (2) the 
stars with middle ages of 10 — 60 Myr (group 2); and (3) the old stars with ages of 60 — 125 
Myr (group 3). 

Allowance for the Galactic Rotation 

We took into account the general Galactic rotation using the parameters found previously 
(Bobylev 2004): Uo = —28.0 ± 0.6 km s“^ kpc“\ a;' = -1-4.17 ± 0.14 km s“^ kpc“^, and 
u'f = —0.81 ± 0.12 km s“^ kpc“^. The Galactocentric distance of the Sun was assumed to 
be Ro = 7.1 kpc, which corresponds to the short distance scale (Dambis et al. 2001). The 
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radial velocities and proper motions of the stars were corrected for peculiar solar motion with 
respect to the LSR using the values obtained by Dehnen and Binney (1998): {uq,Vq,Wq) = 
(10.0,5.3,7.2) km 8"^ 

RESULTS 

A preliminary analysis of the stellar space velocities showed the necessity of setting a 
limit on the total stellar velocity, < 30 km s“\ which we use below. The 

residual velocities U, V, and W were calculated from standard formulas (Kulikovskii 1985). 
We used stars with parallaxes n > 1.5 mas (r < 667 pc). The number of stars used is given 
in the last column of Table 1. Thus, each of groups 1, 2, and 3 always contains a fixed 
number of stars, 114, 342, and 150, respectively. For our sample of stars, the mean error 
in the parallax is 15 — 17%. The mean errors in the proper motion components of the stars 
4.74r/i;cos6 and 4.74r/if,, including the errors in the parallaxes, are ~ 1 km s“b For ~ 40% 
of the stars, there are no data on the errors in the radial velocities (Barbier-Brossat and 
Figon 2000); for stars with available information about the errors, the mean error is ~ 3.5 
km s“^. 

The group 1 stars are distributed exclusively in Galactic quadrant III. The group 2 stars 
fill the solar neighborhood within r ^ 650 pc of the Sun most uniformly. The group 3 stars 
densely fill a compact zone within r ^ 50 pc of the Sun in Galactic quadrant II. We also 
formed a combined group of young (< 60 Myr) stars composed of the group 1 and group 2 
stars. For the latter combined group, the application of all kinematic models proved to be 
possible. 

As the first approximation for the direction of the kinematic center, we take G = 160° 
and Ro = 150 pc (Bobylev 2004). Solving the system of equations (1) for the combined group 
of young stars (group 1+group 2) yielded the following kinematic parameters of the linear 
motion of the stars: {uq, vq, wq) = (0.0 ± 0.7, —12.3 ± 0.7,1.3 ± 0.3) km s“^. Consequently, 
the stars under consideration move with respect to the LSR at a velocity of Vq = 12.3 ± 0.7 
km s“^ in the direction Lq = 270 ± 3° and Bq = 6 ± 1°. Further, 

Uo = —32.4 ± 3.8 km s“^ kpc“^, (12) 

= +93.2 ± 28.3 km s-^ kpc-^, 

= —170.5 ± 107.8 km s“^ kpc“^, 
ko = +27.9 ± 3.8 km s“^ kpc“^, 
k'^ = -122.3 ± 28.3 km s'^ kpc-^, 
k'l = +233.7+ 107.8 km s'^ kpc-G 

lo = 128°, 

(G)d = 2°, 

(G)f = 181°, 

(Ro)d = 180 pc, 

{Ro)f = 410 pc. 

The results of the solution of (12) are presented in Fig. 1, where line 1 indicates the 
derived expansion curve with the estimated parameters fco, and k”] line 2 indicates the 
rotation curve with the estimated parameters cUo, and a;". As we see from the solution 
of (12), the value of Ro = 180 pc calculated using the first derivatives is in satisfactory 
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agreement with the Ro = 150 pc found previously (Bobylev 2004). The expansion curve 
constructed by solving (12) has a velocity maximum of = 4.4 km s“^ at R = 200 pc from 
the center. Based on these data, we determined the time elapsed from the onset of expansion, 
r = 44 Myr, and it is in good agreement with the ages of the stars under consideration. 

Applying the Ogorodnikov-Milne model (Eqs. (11)) to the combined group of young stars 
(group 1+group 2), we found the following components of the tensors M, M+, and M~ (km 
s“^ kpc“^): 


M = 


^ 13.0(3,3) 
-21.5(3,3) 

\ 4 . 6(3,3) 


8.1(i,4) 
2 . 9 ( 1 , 4 ) 
4.3(i,4) 


—4.4(8,6) 
— 7.0(8,6) 
-18.8(8.6) 



( 

13.0(3.3) 

-6.7(i.8) 

0. 1(4.6) 

M+ = 


-6.7(i.8) 

2.9(1.4) 

— 1.4(4.3) 


\ 

0. 1(4.6) 

— 1.4(4.3) 

— 18.8(8.6) 


/ 

0 

14.8(1.8) 

—7.1(5.o) 

M- = 


— 14.8(1.8) 

0 

—4.5(4.6) 


V 

4.5(4.6) 

5.7(4.3) 

0 


(13) 


The Oort parameters are: Aq = = —6.7±1.8 km s“^ kpc“^, Bq = = —14.8±1.8 km 

s“^ kpc“^, Cc = 0.5(M]i[ — M 22 ) = 5.1 ± 1.8 km s“^ kpc“\ Kq = + M 22 ) = 8.0± 1.8 

km s“^ kpc“h In Eqs. (11), the X axis is directed toward the Galactic center (G = 0°). The 
deformation tensor in the principal axes Mq is (km s“^ kpc“^): 


M+ = 


/16.4 

0 

V 0 


0 

-0.4 

0 


0 

0 

-18.9 
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the directions of the principal axes are: 


Li = 153.5 ± 0 . 2 °, =- 1.1 ± 0 . 0 °, ( 14 ) 

±2 = 243 ± 10 °, ^2 = 4 ± 8 °, 

±3 = 81 ± 10 °, 53 = 86 ± 11 °. 


The vertex deviation in the XY plane that we determined from the formula tan2/xv = 
—C/A is IxY = 19 ± 6° and indicates one of the directions of the rotation center; the other 
direction is I = 109° and is close to the direction of the kinematic center lo = 128° that we 
found using Eqs. (1). 

An analysis of the deformation tensor (19), (20) leads us to the important conclusion 
(which cannot be drawn from model (1)) that expansion takes place exclusively in the 
Galactic XY plane. 

The results of applying the statistical method are given in Table 2. As we see from this 
table, the direction of the rotation center for the combined group of stars (Group l±Group 2) 
is determined accurately, /i = 132.7 ± 0.4°. The first axis of the residual velocity ellipsoid 
(the vertex deviation) changes its direction with the age of the stars under consideration. 
This is consistent with the dynamic model of Olano (2001). As Olano showed, the direction 
of the vertex rotation depends on the critical stellar density in the solar neighborhood p/p*, 
where p is given in Mq pc“^. Eollowing Olano’s estimates (p/p* = 1.5), we assume that 
the change occurs counterclockwise in the coordinate system under consideration. Thus, the 
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vertex deviation changes from li = 36° for the sample of old stars (Group 3) to /i = 131° for 
the sample of middle-aged stars (Group 2 ) and to /i = 172° for the sample of the youngest 
stars (Group 1 ). Associating this direction with the direction of the kinematic center, we can 
see that the center does not remain constant. Therefore, its independent determination is 
required for each sample of stars. There is good agreement of the orientation of the residual 
velocity tensor for the Group 1 and Group 2 stars that we found (Table 2 ) with the geometric 
characteristics of the Gould Belt (Torra et al. 2000). For example, for the group-2 stars, the 
directions h = 224 ± 18° and 62 = —22 ± 5° almost match the characteristics of the Gould 
Belt VIg = 275 — 295° and ic = 16 — 20° (here, Qg is the direction of the line of nodes of the 
Gould Belt and zg is the inclination of the disk to the Galactic plane) found by Torra et al. 
( 2000 ). 

For the combined group of stars (Group l-|-Group 2), we obtained the solutions of Eqs. 
(7) in which we used lo = 128° and Ro = 150 pc. Equations (7) were solved for two cases, 
with nine and seven unknowns in the first and second cases, by assuming that the second 
derivatives u'J and /c" are equal to zero. The results are presented in Table 3 and Eig. 2. In 
Eig. 2, curves 1 and 2 were constructed using the data in the upper part of Table. 3; curves 5 
and 6 were constructed using the data in the lower part of Table 3 (seven unknowns). Gurve 
4 in Eig. 2 represents the results of applying the Oort-Lindblad model, which for lo = 128° 
are the following: Ag = 6.5 ± 1.8 km s“^ kpc“^ and Bg = —14.8 ± 1.8 km s“^ kpc“^ (these 
values were obtained from the equations of the Ogorodnikov- Milne model, but the curve is 
essentially the Oort-Lindblad approximation). The time elapsed from the onset of expansion 
for curves 2 and 6 is about r = 120 Myr, which significantly exceeds the estimate obtained 
by solving ( 12 ). 

The results of solving Eqs. (7) for group 3 are given in the lower part of Table 3. The 
parameters derived for group 3 differ significantly from those of young stars: the rotation 
and the expansion have opposite signs. We associate the kinematic peculiarities of the Local 
system of stars that were found by Tsvetkov (1995a, 1995b, 1997) by studying type-A and 
-E stars with these stars. It is the influence of the Group 3 stars that led to the paradox 
that younger OB stars showed a slower rotation than older stars (Bobylev 2004). 

Eigures 3 and 4 show the residual velocities 77, V, and W of middle-aged Gould Belt 
stars in projection onto the Galactic XY and XZ planes. These velocities are residual in 
every sense, because they were corrected from both the peculiar solar motion and the total 
velocity Vg that we found. Thus, we relate the velocities to a coordinate system whose center 
is the center of mass of the Gould Belt that moves with respect to the LSR. 

Eigure 5 shows the residual velocities of stars in the young f3 Pic and Tuc/HorA clusters 
together with TWA cluster stars. 

DISGUSSION 

The kinematic parameters that we derived using the linear Ogorodnikov-Milne model are 
in excellent agreement with the kinematic model of the Gould Belt suggested by Lindblad 
(2000), which is based on an analysis of the results by Gomeron (1999) and Torra et al. 
(1997) for the age interval of Gould Belt stars 20 — 40 Myr. Lindblad’s model assumes the 
proper rotation of the Gold Belt stars with an angular velocity oi ojg = B — = —24 

km s“^ kpc“^ (Aoj = 6.4 km s“^ kpc“^ and B = —17.4 km s“^ kpc“^) and the expansion 
with K = 11.3 km s“^ kpc“^ in the direction of the center lo = 127°. An analysis of the 
curves in Eig. 2 indicates that the linear model faithfully describes the observed velocity 
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field of the stars under consideration no farther than 200 — 300 pc from the kinematic center. 
Our rotation parameters for the Gould Belt are consistent with the model by Olano (2001) 
and with the evolntion model of the yonngest moving clnsters by Asiain et al. (1999). The 
qnestion of whether the f3 Pic, Tnc/HorA, and TWA clnsters belong to the Gonld Belt 
structure is currently being debated. Makarov and Fabricius (2001) directly associate these 
clusters with the Gould Belt; Song et al. (2003) cast doubt on this association. In our opinion, 
a comparison of the stellar residual velocities for all of the clusters {{3 Pic, Tuc/HorA, and 
TWA) in Fig. 5 with onr rotation and expansion cnrves snggests that the clnster stars belong 
to the Gonld Belt strnctnre. The last row in Table 2 gives onr calcnlated parameters of the 
residnal velocity ellipsoid for all 33 stars of the f3 Pic, Tnc/HorA, and TWA clnsters. An 
analysis of /i, 6i, I 2 , & 2 , h, h and a comparison of Figs. 4 and 5 indicate that these stars 
belong to the disk of the Gould Belt. Ortega et al. (2002) performed model calculations of 
the orbits of f3 Pic clnster members in the Galactic field of attraction for an interval of 11.5 
Myr and nsed them to determine the birthplace of the clnster. As follows from an analysis 
of onr rotation cnrve for the Gonld Belt (line 1 in Fig. 2), the Gonld Belt has a snbstantial 
mass. Therefore, in our view, to determine the birthplaces of the f3 Pic, Tuc/HorA, and 
TWA clusters, the gravitational potential of the Gould Belt (^ 2 x IO^Mq) or the entire 
Local system of stars (^ 2 x IO^Mq) mnst be taken into acconnt. 

CONGLUSIONS 

We have analyzed the residnal velocity field of nearby and yonng Gonld Belt stars by 
assnming the existence of a single kinematic center (the center of proper rotation and ex¬ 
pansion). The direction of the kinematic center of the Gonld Belt was fonnd to be: G = 128° 
and Ro = 150 pc. 

Using the Ogorodnikov-Milne model, we showed that the vectors of the radial component 
of the residual velocity field (a positive iF-effect) lie in the Galactic XY plane unrelated to 
the plane of symmetry of the Gonld Belt disk. This led ns to conclnde that the effect in 
qnestion resnlts from the impact of a spiral density wave propagating precisely throngh the 
Galactic disk on the cloud of gas out of which the stars under consideration were subsequently 
formed. In contrast, the effect of proper rotation is closely related to the plane of symmetry 
of the Gould Belt disk. 

We have considered an approach that, in onr view, allows the observed residnal ve¬ 
locity field of the Gonld Belt stars to be described more faithfnlly than does the linear 
(Oort-Lindblad) model. The linear velocities that we found for a sample of stars younger 
than 60 Myr reach their maximum at a distance of ~300 pc from the kinematic center and 
are —6 km s“^ for the rotation and +4 km s“^ for the expansion. Our constructed proper 
rotation cnrve for the sample of yonng stars in the Gonld Belt snggests that its mass is 
snbstantial. This conclnsion is in agreement with the spatial distribntion of faint stars with 
X-ray emission (Gnillont et al. 1998). We argne that the yonng f3 Pic, Tnc/HorA, and TWA 
clusters belong to the Gould Belt structure. 
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Table I. Data on nearby clusters and OB associations 


Cluster/association name 

Age, iMyr 

Reference 

N, 

Group 1 




Ori OB la 

11.4 ± 1.9 

Brown el al. (1994) 

29 

Ori OB lb 

1.7± 1.1 

Brown et al. (1994) 

20 

Ori OBlc 

4.C±2 

Brown et al. (1994) 

25 

Ori OB Id 

<1 

Brown et al. (1994) 

3 

Col 121 

5 

de Zeeuw et al. (1999) 

20 

T\V Hydrae(TVVA) 

< 10 

Song el al. (2003)(f) 

17 

Group 2 




Per OB2 

2-15 

Jimenez and Clausen (1994) 

9 

NGC 2232 

30-60 

Torra el al. (2000) 

4 

NGC 2451 

50 

Carrier et a/. ( 1 999) (c) 

5 

NGC 2516 

30-60 

Torra et al. (2000) 

6 

NGC 2547 

.30-40 

JeRries and Tolley (1998) 

5 

Vel OB2 

20 

de Zeeuw et al. (1999) 

29 

IC 2391 

35-55 

Montes el al. (2(X)1 )(d) 

10 

Trio 

15-30 

de Zeeuw et al. (1999) (c) 

6 

a Car 

<60 

Platais et al. (1998) (c) 

4 

IC 2602 

30 

Luhman (2001 )(c, d, g) 

8 

LCC 

16-20 

Sartor! et al. (2003) (a) 

4C 

UCL 

16-20 

Sartor! et al. (2003)(b) 

67 

US 

8-10 

Sartor! et al. (2003) 

44 

Cep OB2 

.5-10 

de Zeeuw et al. (1999) 

16 

Lac OBI 

15-25 

de Zeeuw e( al. (1999) 

29 

0 Pic 

<30 

Songet al. (2003) 

8 

Tucana/HorA 

<30 

Songt’l al. (2003)(h) 

8 

MaKapoB(XY) 

<60 

Makarov (2003) 

24 

Wiclimann-l 

<14 

Wichmann el al. (2000) 

2 

\Vicliniann-2 (ZAMS) 

« 90 

Wichmann et al. (2003) 

2 

Group 3 




Cep OB6 

« .50 

de Zeeuw el al. (1999) (e) 

10 

a Per(PerOB3) 

« 50 

de Zeeuw et al. (1999) 

46 

Cas Tau 

« 50 

de Zeeuw £>/ al. (1999)(e) 

56 

Pleiades 

70-125 

Soderblom et al. (1993) 

14 

Wichmann (PL) 

Ri 100 

Wichmann et al. (2003) 

30 


Note: (a) 17—20 A\yr(Mamajek et al. 2002): (b) 15—22 Myr(Mamajek el al. 2002); (c) 30—60 Myr(Torra el al. 2000); (d)30 Myr 
(Staufferera/. 1997);(e)50—60Myr(Asiainct al. 1999);(f) 10Myr(Luhman 2001 );(g)30Myr(Randich etal. 1997);(h)30+5° Myr 
(Luhman 2001). 
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Table 2. Principal semiaxes of the residual velocity ellipsoid £ti, 2,3 and their directions /i,2,3, ^1,2,3 


Group 

Age, iMyr 

(Ti, <t 2. 0-3, km s 

h.h.h 

5i, 52,53 

1+2 

<60 

9.1 ± 1.9 

132.7° + 0.4° 

-6.0° +0.0° 



6.8+ 1.9 

225° + 27° 

-19° +4° 



2.8+ 1.9 

206° + 27° 

70° + 3° 

1 

<10 

7.6 + 2.3 

172° + 1° 

-10° + 2° 



7.2+ 2.3 

263° + 42° 

-2° + 6° 



3.5 + 2.2 

186° + 42° 

80° + 6° 

2 

10-60 

9.4+ 2.0 

131°+ 1° 

-7° + 1° 



5.5+ 2.0 

224° + 18° 

-22° + 5° 



2.7+ 1.9 

20.5° + 18° 

66° + 4° 

3 

50-125 

10.3 + 2.3 

36° + 1° 

4° +0° 



6.1 +2.3 

127° + 10° 

17° + 4° 



2.5+ 2.3 

293° + 10° 

73° + 5° 

R Pic + Tuc + TWA 


5.6+ 1.7 

47° 

34° 



3.5+ 1.7 

126° 

-16° 



2.1 + 1.6 

194° 

51° 


Table 3. Kinematic parameters of the Gould Belt stars calculated from Eqs.(7)with the assumed = 128°, Ro = 150 pc 
for group 1 +2 and lo = 36°, = 150 pc for group 3 


Group 

Age, Myr 

UG, CQ. 

Pg, Tjg, Rq 


fro, K. fr" 

1+2 

<60 

1.1 + 0.4 

11.9 + 0.4 

-35.3 + 3.6 

12.1 + 2.7 



-11.8 + 0.4 

275° + 2° 

132.6+ 25.7 

-9.0 + 13.2 



1.3+ 0.3 

6° +1° 

-292.5 + 87.9 

-47.6 + 29.9 

1+2 

<60 

2.3 + 0.3 

12.0 + 0.4 

-23.1 + 2.2 

14.0 + 2.2 



-11.1+0.3 

285° + 2° 

31.3 + 6.5 

-27.3 + 0.5 



1.3+ 0.3 

6° +2° 

- 

- 

3 

50-125 

0.7+ 0.8 

11.2 + 0.8 

.54.7 + 9.2 

-23.2 + 9.2 



-11.2 + 0.8 

273° + 4° 

-103.5 + 42.6 

123.4 + 42.7 



-0.1 +0.5 

-1° +2° 

- 

- 


Note: The linear velocities uq. vq, wq, and V'g are in km s the angular velocities iCo and fro are in km s ' kpc ^ the derivatives oco 
and fr-' are in km kpc~^; and the derivatives u,'" and fr-" are in km kpc“®. 
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Fig. 2. Expansion (2, 6) and rotation (/, 4.5) curx'es of Gould Belt stars with ages <60 .Myr(group 1 + group 2); the position 
of the Sun is marked (<?). The parameters of the cun-es were obtained for the centerw'ith lo = 128° and Ro = 150 pc. 
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Fig. 3. Residual velocities U and V' of Gould Belt stars with ages <60 Myr( Group 1 + Group 2) in projection onto the Galactic 
XY' plane. The Sun is located at the coordinate origin. The circle marks the kinematic center with= 128° and Ro = 150 pc. 



X, pc 

Fig. 4. Residual velocities U and U' of Gould Belt stars with ages <60 Myr( Groups 1 and 2) in projection onto the Galactic 
XZ plane. The Sun is located at the coordinate origin. 
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npc 

Fig. .5. Residual velocities of () Pic and Tuc/HorA cluster 
stars (the positions of the stars from both groups are 
indicated by the circles) together with TWA cluster stars: 
U and V' in projection onto the Galactic A')' plane (a); U 
and If in projection onto the Galactic XZ plane (b); V 
and IF in projection onto the Galactic YZ plane (c). 
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